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SITIMARY 

Tlje linear trnnHonlc perturbation integral equation previouBly derived 
for nonlifting airfoils Is formulated for lifting cases. In order to treat 
shock wave motions, a strained coordinate system is used in which the shoetc 
location is invariant. The ;angencv boundary conuitions are either formulated 
using the thin airfoil approximation or by using the analytic continuatieM 
concept. A direct numerical solution to this equation Is derived in contrast 
to the iterative scheme Initially used, and results of both lifting and 
nunlifting cx.imples indicate that the method la satisfactory. 


I. INTRtn»urTU»N 


AN IMt’i'RTANT piitbtom In iUM»ulvn.imli*» Im th«* .ici'nr.ito t Inn nl 

tiu' pr**Hauro» on an alrfntl that Ik nHt'illalInK hai innnii'ally wtiJt 
Miiutll amplltiulo In a tranMonlr In tIu* limit rnrn I rn- 

qnonry thlH proMom rt*klui'OH to that of an airloll that hait utuli'i'** 
goiu* a perturbation equal to the magnitiule ot tlu> OHei I lat ionM. 

Thitt Hlmpler Hteadv perturbation in theretore a gooil starting 
point to ilerivlng a fundaroiMital approaelt for unsteailv flowH or to 
gain experlonee in »leveloplng num«.*rteal proeeilureti . 

In tranrtonie I lows where there Im a kliiieont inui t v , eortaln 
problems ean arise in a perturbation soliit ion in the region bouiuleii 
by the extremities of the slu’ek nH>vement inilueeii bv the perturba- 
tion. I’erturbatlon solutituis assuim' that elianges in (tie t low 
vari.ibles are of the same otiler of imignitUiie as the lirivlng pertur- 
bation. In thi‘ region traverseil bv t lie shock tiurlng tlie perturba- 
tion the flow variables can jump t rum ,i pre-shock value to a post- 
shock value, aiul henct* c.innot be consiiiereki siwill per t urbat ik'tis . 

Although little attention h.is been paid to the ste.ulv pertur- 
bation pik'blem finite klltfereiui' Sk'Iut Ituis are available for the 
simll.ir hainxMilc.il ly ilecomp»>sed unsieadv tr.insonlc t Kkw pu’blem. 
I’reseut ly .ivallable metluvls i*f solving the h.irmonical Iv decompk'seil 
problem Incliule the ilnlte illflerence calt'oial iiui ot Ir.icl et al.* 


In thoai* iMHltoila no account la t.ik(*n of 


2 

aiul Woutliorhl 1 1 ot nl, 
tin* aliock naUion, the aliock bolnit aitaumoii to rom.iln unmciU ia 1 1 y at 
Ita steady state location. A serious dr.iwback ot thusi procedures 
Is the occurrence of a severe numerical Instability In tho relaxa- 
tion procedure used to solve the dilterence uquatlons. This in- 
stability occurs at a critical Mach number depondent frequency bevond 
which the method diverges. A means ot counturlnp the difficulty 
Is (t^ven by Hafez et al.;^ however, the rate of convergence is slow. 
An alternative nsMns of avoiding the Instability Is by direct 
solution of the difference equations, l.e., without Iteration. 

However, the large number of mesh points makes this procedure not 
feasible especially for practical three-dimensional flows. All ot 
these methods are shock capturing algorithms. Shock fitting pro- 
cedures are given In Kef. 3 for the harmonic decomposlt Ion. A 
method of directly Integrating the unsteady transonic sm.il I dlstur- 
bance equation Is the ADI method of Hallhaus and Coorjian . 

An alternative to the finite difference solutions for transonic 
flow Is the Integral equation method^’^. The basic advantage ot 
the integral equation method Is that considerably fewer points are 
required for the discretization. In the present context this implies 
that direct inversion of a discretized equation Is a more feasible 
proposition than for the finite difference methods. Also from 
evidence of the steady version of the Integral equation diethod an 
Integral equation solution of the perturbation or unsteady problem 
would he computationally faster than the direct Integration method of 
Kef. 4. Furthermore the Integral equation method gives automatic 
shock fitting. The main difficulty cf the perturbation problem Is In 
the treatment of the moving shock location. 


- 1 - 





A irtf.iiiii of oviMcomIng iIiIn illfllcultv I* nivon bv Nixon In 
whli'h a Ntralni'il ooiTillnato Mvatom la iiaoil. In tboao now i'ooi'iIIiuKom 
tbo abook Inoatlon ronu’ Ina conatant ami tho maitnltmlo i>( (lio 
alralnlnit, whloh la rolato«I to tba abook movom»*nt , la loun*! aa 
part of tbo aolntlon. In roforonoo (7) tbo roaultlntt linoar pi*r- 
turbatlon oijnat ton for nonlltttnit flowa la aolvoil natiiK tbo tranaonlo 
tntoitral oquat ton (S) molluul wltb an lloratlvo pn*oo»biro. Tbo 
nkixnltUilo of tbo abook nx*vonH’nt , ami bonoo tbo illatortlon of tlio 
ooorkllnattf ayatom, la founil by onlorolng tbo rogularlty condition 
wboro tbo flow la ounttnuoua ovorvwboro apart from tbo abook wavo. 

Slnoo tbo baalc oquatlon la Itnoar It aoeraa that a moro dlroot 
aolntlon of tbo portnrbatlon onnatlon la poaatbto, and It la tbla 
approaob that la Invoatlxatod boro wltb tbo prinotpal objoottvo 
bolnx a tuturo oxtonalon to oaolllatory dlaoont Imioua tranaonlo 
flowa. Alao, tbo portnrbatlon Intogral oqnat ion for Itfttnx flow 
la dorivod. 

Tbo Intoxral ounatlon la rodnood by qnadratnro to a avatom of 
1 Inoar algobralo oqnatlona wbiob Inoorporato tbw roxnlarlty oondt- 
llon of oontlnnona flow oxoopt at tbo abook. I'boao o()natfona aro 
aolvod by t'.anaalan ollralnatlon. Satlafaotorv roanlti for llftlnx 
and nonllftlnx portnrbat Iona Indicate that tbo prooodnro o.in bo .ippltod 
wltb oonfldonoo to tbo mtiob moro Important oaolllatory probloma. 

Snob an application la reported In roforonoo fS). 

rbo "baao" aolntli>n naod In tbo oalonlattona roqnlroa a noinuil 
d laoont Innona n'proaontat ton of tbo alu'ok. In tbo proaont work, 
tbo baao aolnt Ion la oompntod by a abook oaptnrlnx finite dllfoionoo 




program with a fine meah since this is the miiat convenient meana 
available. The nornuil, diacont inuouH shock requirement 1 h obtained 
by suitable entrapolat ion of the velocities outside the shock 
capture region. 


2. BASIC EQUATIONS 

In a Cartesian coordinate system (x,i) with the origin at the 
leading edge of the airfoil and the x-axls aligned with the air- 
foil chord, the transonic smal 1-disturbance equation for a free- 
stream Mach number can be written as 


*xx ^zt " 2 ^^x^^x 

where, if 6 ■ (1 - y 1" the ratio of specific heats, 

and k ■ kCviMo,) is a transonic similarity parameter, then ^ is 
related to tlie perturbation velocity potential i by 


(>(x,z) ■ i(E,*) , 


and the coordinate system (x,z) is related to (x,z) by 

X ■ X 1 


( 2 ) 


z ■ 0z I . 

In this transformed space, the perturbation velocities (u,w) in the 
X- and z-directlons, respectively, are related to the physical 
perturbation velocities (u,w) by 


u 


w 


li 

9x 


21 

dz 



k 11 ’ 

► 

A. 

6^ 3z / . 


(3) 
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TW K'uiuUiry coiuiltionM arc that 

a) the* fli^w at the airfoil Hurtai'o ia tanttontial ti« the 
airfoil Hurfaee, 

h) the veloeitv ia finite an infinite tiiatame from the 
airfoil, anil 

c) that for a Hubaonic trailinK eilge the Kiitta comlilton of 
finite preaiuire at the trailing eilgr miiat be Hattafteit. 

If the equationa of the upper and lower aurfacea of the airfoil 
are niven bv i •• lu(S) and X ■ X|,(x)t reapect ively , then the 
tanttenov boundary condition in the tranaform«*d coordinate avatem 
(x.s) can be written aa 

wix.jtj) • Z|^(XI^I ♦ ^ u(x,*^)J - \ 

where if a ia the angle of attack then 
T k 

V • — J vl 

- k 

r^ix) • t<Xj,ix) , Z„(x) - 

Xj (x) • ^X( (x) , ^ *1 

For thick airfoila, the tangenev boumlary condition on the 
airfoil aurface, equation (4>, can be replaced bv an equivalent 
boundarv condition on the airfoil chord line, x • *0, bv using the 
idea of an analytic cent inuat ion (V) of the external flow inside 
the airtoil. This is obt.iined by using a Tavlor series to expand 




-n- 


w(<, I ) and w(x,z, ) tu th« chord lln«. Thus, on uslnK squatlon (I), 

il L« 

equation (4) gives, to second order in magiltude of the thickness/ 
chord ratio of the airfoil. 


- -A + Z (x) 

1 + u(x,* ) 

* (x) 

u(x,* ) 


u 

X 

k u 

u 

u 

1 


I 2, . 

- J u (X.I^,) 


w(x,-0) • -A + (x) 

X 


(6) 


1 + ^ u(x,*j^) 


•* «t‘»> 


u(x,Xj^) 


-iu^x,*L)J^. J 

The usual thin airfoil boundary conditions can be recovered from 
equation (b) by neglecting all the second-order terms. This analytic 
continuation device is necessary if the Integral equation method (S) 
is to be used in its commonly derived form in which the boundary condi- 
tions are required on the chord line rather than on the airfoil surface. 

There is an apparent inconsistency in the use of a second order accu- 
rate boundary condition like Eq. (6) with a first order differential equa- 
tion, Eq. (I). The small disturbance approximation breaks down in the 
region of the leading edge and leads to singular behavior. Generally 
this singular behavior is removed by including certain (fairly arbitrary) 
second-order terms. The boundary condition of Eq. (6) can be regarded as 
simply a more consistent way of introducing the necessary second-order terms. 

The weak shock Jump relations for equation (1) are 




tan 



J 


( 7 ) 
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win* rt* 


* 

[ ] lienoti'H .1 Jump acroHM a aluuk wav»* an»l la tlu* .iiikIv tin* 

alu»ik makra wlih tho a-axla. If tho ahuck wave can be aaaumfJ nornuil to 
the free at ream, then 



f 

■ 0 


IH) 


The problem uiiiler conalderat Ion la the calculation of the effect of 
a aimill change In the boundary lerma wlx,«^>, wlx.Xj > of equation (4) 
on the preaaure dlatr Ibut Ion. Tliua, let equation (b) be written aa 


w(x,>0) - Wq(x,MI) + iW|(x,+0) 
w(x,-0) « w^tx.-O) + iWj(x.-O) 


l‘») 


where C la a anuill paranu'ter. 

The problem la therefore to aolve equation (1) aubject to the 
boundary condition equation I**) aaaumlnR that the aolution of equation 
(1) aubject to the boundary condltlona 


wlx.+O) " Wq(x,-H)) 
wix,-0) • w^lx,-0) 

ia known. 


UO) 


It ia ahown in reference (7) that If ahock waves are present In 


the flow then a distortion of the toordlnate svstem la necessarv In order 
to correctly formulate the perturbation problem. The shock location Is 
invariant In this distorted coordinate system. It Is a basic rest rlit Ion 
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of IIk* theory thjt ahoek wave* r<inuot he Konerate ' tJentroved during 
the perturb.it Ion. It 1« .iluo .l•llu■ed• (or •implicltv, thut all ahock 
w.tve« o.in be eoneldered normal to the free «tie.im. Thin me.ina that only 
the at re.iRiw(iie eoordlnate x need he atralnad. following the general 
line of attaek ot refereiu-w (7) let there he ahoek waves in the upper 

half plane, and shock waves in the lower half plane. 11 x* denotes 
the strained streamwise coordinate, then let 


X • X* ♦ iXj(x') ♦ ••• 


Ul) 


where, following the suggestion of reference (7) 

Xi(x') • 5x^ Xj(x') 0 ^ x' < I 

i-1 ‘ 

Xjix') - 0 x' < 0. x' ' 1 

where N^. - N ^ N, . 

Sul 

In equation il«) 5x^ is the sk'vement of the 

* j 

x^ix') is some straining function. If the velocity 
expanded in the S'^riss 


(12) 


ith shock wave and 
potential 4'(x,i) is 


- ^q(x’,*) ♦ L«>j(x',*) ♦ ••• 


(IJ) 


then on substitution of equations (12) and (13) into equation (1) and 
equating coefficients ot i , the equations for and are obti ined; 
thus 


0 , , 
XX 


tt 


Jo 


X X 


(U) 


with the bounda r y cond i t ion 


(x‘ ,‘0) - w^^(x* ,t0) 

c 


(l^) 



ilb) 
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with Ihv boundary condition 


(x*,*0) • w (x'.tO) ♦ x.(x’)Wg tx',‘0) . 

* '^x* 


(17) 


It lx aaxumod that tho solution to the problem doflnud by equation 
(la), equation (1)) Is known. 

The normal shock Jump relations for equation (lb) are 


♦ l -♦! 

X X 


•k' ■ ' *«.•) 


['>if • 0 


(i«) 


F.q. (16) and Eq. (18) can be written In Intuitral form using Green's 
theorem. The domain of Green's theorem Is the entire flow fluid with the 
exception of the slit z ■ any shock waves In the flow, and the field 
point (x,z). In this procedure the line Integrals around the shack waves 
vanish because of the shock conditions, Eq. (18). If the general Ideas 
of Ref. (6) are followed then ..he Integral equation for the system defined 
by equations (16), (17), (18) Is as follows. 

Uj(x',z)[i - Uq(x',z)] - Ij^(x'.z) -E ljj(x'. ,z,x^) 






(I**) 


where 


1-1 


Ujj(x',z) - (x',z) ; Uj(x',z) • (x'.z) , 

^x' x' 


(20) 


and xA denotes the location of the '"base" shock waves. For z ^ 0 


r* Aw (C)(x' - 4) , f 


, , Aw (C)(x' - O 

j ,7^ ,,2, 2, 

0 l(x - C) -E Z J 


Auj(C)z 




— d^ (21) 


lg(x' ,z,x^) 


-hfl 


(x* .^;z,;)uj(4,:.)uQ(^,;)dC d; 


( 22 ) 


- 1 (x' , z, x'^) 
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/, (*' KiCOwj (0(x* - n 

‘ \ 1(X’ - ry ♦ t^] 


llf 


♦ X< 


, (x' > [2 Uo(k’ .«) - j u^^(x',*)J 






'.C;z.C)^Ug(t.c^ - y U(,?(^c>]xj (f.)jdC d j 


(23) 

wh»*ro ihf oporatur A (h dof Inod for a function f(C,'.) bv 


Af(0 - f(C.+0) - f(C.-O) . 
Kor c ■ *0, In aquation (H) 


(24) 


Ij (x’.tO) 


± ,, f‘ “i‘‘' I f i'\, 

2" J (x* - o ^ X \ X* / J (X* - O \l - r/ 

•'o •' n 


(2S) 


Ic(x’.f 


.^o.xj) . . J 

•'o 

> Is^(x'.O.x') 


It.(C,x^) 

oT^'o 




(2n) 


/ (C.Xj,) . X I 

i,^<x’,.o,x;, - , I at 

•'o 

•*• " I 
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JJ (x’. <.;0. (,)(*!, - y{f.,^)U\r. 


mul 


K ( C • i# ^ 


u,(r., ♦«»)»(, (t.H)) . c ' 0 

» 1 , (^..-O)Up(t.-O) . t ‘ 0 


IrIx ', k ||.> - Wj ( x ’) ♦ jl ^; 

" o ^ 1 \ 

- - y - ((..Oj ‘«t 


2 Ho ( C ...' - y 


wlu> 


v(t.O - 


. ; > 0 




” 2 "0 

Ttu* Inli'Ki'ol ovt'i' S tM (IH liu’il l>v 

/r'"'’ r r 

I If «u » ir . - Urn I I I K lit ac fl | y at at 

if J'.-y- •'« J.'«, ^0 


I I -'I I 

•'x'. • • •'0 » 


F ar. at 




w»un»* Me i»n»l Mi roprfMfnt nil nluu-kN In llu* uppnr rtiul lowr 
^n ‘ I, 

ImiH plam'M I vt»l v, niul 

w^(0 - \ |wj(C.^O> ♦ w,(t. '0)1 

Wj(m') - 2 .♦«’> *■ • 

It can he Hc«>n that tM|uat ion (I**) will nlve an Infinite value 


ol U|(m\m) when UpfM'.c) • I (the aonlc line) iinleNH 


«* 

S 


«Mc 1, (m'.I.Mc) 

”1 ’ i ^ 


M*“M- (*) 

Op 


I l| (M* ,1-) ♦ lp(M’ ,*,Mj.> I 


m'-m« (*) 
I’ 


Ol> 


where m' " M^ (*') (p ■ ItNo) »leiu»tea the locat li»n ol the pth 

Up 

Monlc line. II there are hIiocWh In the flt*w, then equal ii>n ()l) 

HlveH the required equal lonn for the Hlux'k movementH 

i<M„ (1 • l,Nc). 

H, h 


1. Nl'MI-RIl’AI. SOI.UTION I'KOOKIHIKK 
The flrnt «tep In th«' aolul Ion of the Inteitral eipiat Itui (l‘>) 
(h t»> dlacretlje the IntegralM. Tin* technlquea uM»*d to diHcretl*** 
ar< outlined in J . Mrleflv, the flow field Ih divided 

Into Mtrlpa parallel to tlie M-aMlH, a« ahown in Kin* 1 1 -uni the 
variation of the flow variable part ol the Integrand, for example 
Uj(t,CK Im approximated In terraa of valuea on the atrip edjter 
The Held Inteitrala c.in th«>n he integrated with reapect to C to 
give a line Integral in terma of valuea on the atrip edgea, lor 
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i*x;implc< Ug(C,Zj)Uj (C.c^) (J ■ l.N^) wJuto Z| J In* |th 

Hlrlp film* rtiu! Njj 1 h tiu* mimbiT of HtripH. The roNultiit)t line 
liiti'ttral Ik then ovjilu.itoii by qiMilr.it uro in tormH of Hporlfiod 
valuoH of the MlroamwiHo variable wboro Nj^ in tlu* 

mimtuT of X Mtatii'na. The valuea of u^(X|^,Z|), ii|(X|^,Z|) are 
m>w denoted by Ug , «j^(t ■ l,N) reapeot ively , where 

t - k + (J - l)N^ and N - N^N,, 

Uainit this notation the inteKral equ.it ion (I'f) can be written 


aa 


U, (I - Ufl ) • 1 


whe re 


>-T, * ^ "•."•l * £ ■‘•.“".'''j * S 


(32) 


*f.,j • 

Nx 

T ■ *1 “ “l /'I • 

'•^J ^3 7^1 ‘l ‘j 

The matrlcea *'lj* ^ij ‘’valuated uain>; the Ideaa and 

Influence functiona i;lven in , and the ^x^ (i - 1,N„) 

H J W 

are given on eacli atrip edge by equation (31). If a auperacrlpt 
(k) denotea values on the kth atrip edge then from equation ()’) 




) 


( n) 


(k) (k) 

where are the valuea of Ij , I^, respectively, at tlie 

ptli Monlc line on the kth strip edge. The matrix is tlie 

‘ip 
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vrtlf.c of Ij. al tin* ptb Monlc Ilu»’ on the kth Htrip oiIko aiul 
SiihNtUutlon of oqn.it ion (1i) into o(]uat Ion (iJ) nlvoH 



wiioro valiioM at tlu* pth Honli* Itiu* on tlu> kth 

Htrip (o'norally tlu*Ho valiioN aro foiiiul hy int orpolat ion 

from tlu* Htaiuiani valnoH a^|. 

Kquat ton (14) can lu* written in the simple form 


A,tu,^ - H, 


(IS) 


where 



( Ih) 


( 17 ) 
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Ajj is the Kronecker delta, and 



Mj ■ 1 

1 X 

(38) 

Wl • 0 

1 2 . 


Note that the nuitrix Aj^ 

ia singular if the 

x^ coincide 


with a Honic point. Honco the dlatributlon of tho computational 
moah ahould be auch that the aonic pointa on each atrip edge do 
not lie on a meah point. la a large denae matrix and in the 

present scheme the set of equations (3S) ia solved by (•auaaian 
el Imination. 

ikiving obtained the Uj^ the total velocity u.j.(xj,*j^) ia 
given by 

u^fxj.*,^) - Uj,(x[,Xk)jl - ^’‘SpXp^.txp + cu,(xj.z^) (3Q) 

P-1 * J 


wliere 


>*1 


Ns 

i| + e ^ fiXa^ipCxp 

p-1 


(40) 


and x^ ia found from equation (33). 


It can be seen from equation (3d) that the total contribution 


of the perturbation variables ia 

r N. 


CUj(x|,Zj^) - CUj(xpZj^) - 


L p-1 




i)juo(xpr,^) . 


(41) 
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ar«» Holvod 


In the preceding section the ter.M and 
separately and it Is possible that cU|(xpZ|^) is small , wliereas 
its components rU|(X|,r.|^) and 


- Ns 

_ P-> ‘ 


may be relatively large. This situation can lead to large errors 
in cU|(x|,Z|^) even though the relative errors in its component 
parts are snuill. Consequently, it has been found in some cases 
that an alternative formulation in which rather than 

Uj(xJ,Z|^) is the dependent variable is more accurate. Thus 
equation (32) is replaced by 


where 


N 


Apart from the definition of If the numerical treatmeiu is 

PJ 

identical to the previous scheme. Ihiving found u, the total 

‘i 

velocity is given by 

u^Cxj.r.,^) - Uq(xJ.Zj^) + 

and xj is given by equation (40). 
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A. RESULTS 


Two exanpIcB of ptTturbat ion fliWH rniculatod by llu* pruHcnt 
method are shown In Figs. 2 and 3. The Input for these examples, 
the base solution, is calculated by shock-capturing finite differ- 
ence method using the transonic small disturbance equation, with a 
fine computational grid being used. The necessary nornuil shock 
solution is then estimated by the use of a suitable extrapolation. 
In both examples the transonic parameter k is given bv 

k - (y + 1)M„2 . 

For a single shock wave in each half plane the straining function 


il(x') 


■ '‘s> 


is used. 

The boundary conditions used in the first example are the 
usual thin airfoil boundary conditions. In the second example, two 
methods of eliminating the leading edge singularity are used; 
namely, the use of the analytic continuation boundary conditions 
of equation (6) and the use of the well known Riegels factor. In 
the latter, all the terms in equation (26) and equation (27) 
involving the factor [(1 - x)/x]*^^ are divided by 



which removes the x”*/^ singularity for round-nosed airfoils. 

In Fig. 2 the perturbed ^low around a lOX biconvex airfoil at 
zero angle of attack for two different free-stream Mach numbers Is 
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HluHi'n. The b4Mc Mach number ia O.MOM niui tn the preHenl eem|Hil<i> 
t Iona thin airfoil bouiulary eotulitlt*nH are uHed. KlgureM .a and ?b 
ahi*w the preaMure dlatrlhut ion for • 0.81H and • 0.8J8, 
reapeet Ively . Alao alu<wn are the reaulta ot dlreel ealeulat Iona. 

It can be aeen that the preaent mt'thod ittvea kihhI agreement for 
the loM*r Haeh number but that there la poorer agreenu*nt at the 
higher Hiioh number. Thla dlaerepanev la attributed to the linear- 
Ixatlon proeeaa leading to equat Ion (Ib> being no longer a valid 
approxlMbit Ion. 

in Pig. 3 the flow around the upper aurfaee of a NAl*A^•4A00^ 
airfoil at M.,. ■ 0.8S4 at two different anglea of attack la ah«>wn. 
The baae fl»»w la nonlifting. In Pig. la the preaaure dlatrlbutlon 
at 1/4* angle ol attack la compared to tlie reaulta of a direct 
finite difference calculation. It can be aeen that the "analytic 
cont Inuat Ion" reaulta give a better agreenuMit with the direct cal- 
culation for the preaaure than the "Kiegela factor" calculavion. 
However, the ahock location predicted bv the "Kiegela factor" method 
la In better agreement than the "analytic cont liuiat ion" method. 

The probable reason la that In the region of the ahock nu«t ion the 
"Kiegela factor" method givea eaaentlallv the aame linear thlt< 
airfoil boundary conditlona aa for the direct-calculation, whereas 
tl»e analytic continuation boundary conditlona are nonlinear. It 
can alao be aeen that In both examples there la a discrepancy 
between the present reaulta and the direct reaulta especially In 
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thtf leudlnK-i*(lK«^ rcRion. Thin Im probably due to the different 
mi'thodH of treating the leading-edge region in tl«e prenent method 
(analytic continuation or Ricgelit factor) and in the direct aolu- 
tion (thin airfoil boundary condltlona with meah adjuatment). 
Similarly, the rcaulta of the higher angle of attack, a ■ 1/2*, 
shown in Fig. 3b, are not in auch good agreement aa regarda ahock 
location. Again, thla la attributed to the linearization proceaa 
being no longer a valid approximation. 

Generally therefore, the preaent method will adequately com- 
pute the perturbation flow for both lifting and nonlifting caaea. 
For lifting flowR aome dlaagreement with direct rcaulta may be 
expected becauac of dlfferencea In the treatment of the boundary 
condltlona. Alao It would appear that for ahock movementa of more 
than 5X-6% of chord the linearization of the equatlona la no longer 
a valid approximation. 


5. CONCLUSIONS 

A direct numerical technique baa been devlaed for the aolutlon 
of perturbations of discontinuous transonic flows for both lifting 
and nonllftlng configurations. Results of two cx.'imples indicate 
that the method is satisfactory provided the perturbations are 
small. The main purpose of the present paper is to test the 
numerical approach for future use in unstef.dy transonic flow 
problems. 


- 20 - 


Referrnceii 


1. Traci, R. M., Albano, E. i)., and Farr, T. L. , "Small DiaCurbance 
Tranaonlc FIowh About Oaclllatlng Airfoils anu Planar UingH," 
AFFDL-TR-75-10, 1972. 

2. Weutherhill, U. H., Ehlers, F. E., and Sebastian, J. H., "Computation 
of the Transonic Perturbation Flow Fields Around Two and Three 
Dimensional Oscillating Wings," NASA CR 2S99, 1975. 

3. Hafez, M. M., Rlzk, M. H., and Murman, E. M., "Numerical Solution 
of the Unsteady Transonic Small Disturbance Equations," ACARD 
Specialist's Meeting on Unsteady Airloads in Separated and Unsteady 
Flow, Lisbon, 1977. 

4. Ballhaus, W. F. and Goorjian, P. M., "Implicit Finite Difference 
Computations of Unsteady Transonic Flow About Airfoils Including the 
Treatment of Irregular Shock Wave Motions," AIAA Paper No. 77-205, 

1977. 

5. Nixon, D., "Calculation of Transonic Flows Using an Extended Integral 
Equations Method," AIAA J 15 , 3, 1977. 

6. Nixon, D., "Calculation of Transonic Flow Using Integral Equation 
Methods," Ph.D. Thesis, University of London, 1976. 

7. Nixon, D., "Perturbation of a Discontinuous Transonic Flow," AIAA J 10 , 
1, 1978. 

8. Nixon, D., "Calculation of Unsteady Transonic Flow Using the Integral 
Equation Method, AIAA Paper 78-13, 1978. 

9. Nixon, D., "An Alternative Treatment of the Boundary Condtlons for the 
Flow Over Thick Wings," Aeron. Quart. 28 , Part 2, 1977. 


21 



Figure 1.- Arrangement of strips. 
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DIRECT CALCULATION 

— PERTURBATION THEORY 

(ANALYTIC CONTINUATION! 

— O— PERTURBATION THEORY 
(RIEOELS FACTOR! 

— O- BASE SOLUTION 
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FlKurt* 3.- I’r»*s8urc distribution around a NACA (i4A00(> airloil 
at M» ■ 0.854; baso solution a ■ (!® , 


